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Disclination Unbinding Transition in Quantum Hall Liquid Crystals 
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We derive the the long-wavelength elastic theory for the quantum Hall smectic state starting from 
the Hartree-Fock approximation. Dislocations in this state lead to an effective nematic model for 
T > 0, which undergoes a disclination unbinding transition from a phase with algebraic orientational 
order into an isotropic phase. We obtain transition temperatures which are in qualitative agreement 
with recent experiments which have observed large anisotropies of the longitudinal resistivities in 
half-filled Landau levels, lending credence to the liquid crystal interpretation of experiments. 
PACS: 73.40. Hm, 73.50.Jt, 73.20.Mf, 64.70.Md 



Recent experiments M~m m high mobility two- 
dimensional electron systems (2DES) have revealed re- 
markable new phenomena in the transitional regions be- 
tween the different plateau of the Hall conductance. In 
particular, striking anisotropies and non-linearities in the 
magnetotransport were observed for Landau level (LL) 
filling factors near v = n + 1/2, for n > 4, corre- 
sponding to partially filled LL indices L > 2. This 
anisotropy tends to align with the crystalline axes of 
the sample, but can be reoriented by the application 
of in-plane magnetic fields y,|| , and resistance ratios as 
high as Rxx/Ryy ~ 3500 have been observed ||. This 
anisotropic behavior has been attributed to the forma- 
tion of a striped phase. A unidirectional charge density 
wave (UCDW) had been predicted several years ago J7|] 
for nearly half filled high LLs; exact diagonalizations for 
systems of up to 12 electrons B corroborate this picture 
for L > 2, and many experimental results can be qualita- 
tively understood under the assumption of a UCDW. The 
presence of stripes has already been directly observed in 
a large class of low-dimensional, strongly correlated elec- 
tronic systems , and the present experimental evidence 
in quantum Hall devices is compelling, even if still some- 
what circumstantial jlOf. 

Due to the similarities of the UCDW state with a clas- 
sical smectic liquid crystal, these states have been dubbed 
quantum Hall smectics by Fradkin and Kivelson pLH . 
In two dimensions thermal fluctuations destroy the po- 
sitional order jl3|, but the system should still exhibit 
anisotropic transport as long as there is some remnant of 
orientational order (algebraic order in the quantum Hall 
nematic) E3]. As the temperature is increased, the alge- 
braic orientational order will disappear in a Kosterlitz- 
Thouless (KT) disclination-unbinding transition fll5[ . 

To study this process we have mapped the interact- 
ing electron system (in the Hartree-Fock approximation) 
onto a classical smectic (the UCDW). We then consider 
the role of thermal fluctuations (phonons and disloca- 
tions) in reducing the order from smectic to nematic at 
larger distances. Without the use of any fitting parame- 
ters, and using only experimentally accessible values for 
the electron density and the width of the 2DES, we are 
able to estimate values for the disclination unbinding 
transition temperature, which are in qualitative agree- 
ment with the transport measurements. 



(I) Hartree-Fock approximation for the charge-density- 
wave state. — In order to study the energetics of a charge 
density wave (CDW) in the 2DES we closely follow the 
strategy developed in Refs. Pq-Ol, and use the Hartree- 
Fock (HF) approximation, which corresponds to the as- 
sumption that the electronic state can be described as a 
Slater determinant of single-electron states. In the Lan- 
dau gauge, A(r) = (0, Bx,0), and the eigenstates of the 
non-interacting problem are 
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where a, a, n and xq indicate the electric sub-band index 
(due to the confinement in the z direction), spin index, 
LL index, and guiding center respectively; lb = (h/eB) 1 / 2 
is the magnetic length, L y is the length of the system in 
the y direction, and H n are Hermite polynomials. 

Since the electric sub-band splitting is very large 
(about 9.8 meV in the sample of Ref. Q), in what follows 
we consider only states with a = 0. The Coulomb inter- 
action between the basis states above can be replaced by 
the effective interaction [n6|- 
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where k is the dielectric constant of the semiconductor 
(~13 in GaAs/AlGaAs), with the matrix element 
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which may be expressed in terms of associated Laguerre 
polynomials |16|-[18|. Since the anisotropic states occur 
for moderately weak magnetic fields, the effect of a CDW 
on the valence LL is to polarize the fully occupied LLs 
below. This polarization may be accounted for with an 
effective dielectric constant e(q), which can be calculated 
in the random phase approximation (RPA) |17|-|l9| . This 
effective interaction greatly simplifies the calculation, as 
we only need to consider states within the valence LL for 
the determination of CDW energies. 

In the absence of LL mixing, the state of the system is 
uniquely specified by the particle density function plpOf . 
The energy per electron in a CDW state at a fractional 
filling v* is given by plj 
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FIG. 1. Dependence of the average energy per elec- 



tron state E uc 

(n e = 2.67 x 10 11 cm -2 , z rms 



for various filling factors of Ref. 
58.3 A). 
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where A(Gj) is the Fourier coefficient of the occupation 
number at the reciprocal lattice vector Gj and the ker- 
X(q) with the direct and exchange 
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In the UCDW state, we have Gj = e^ G\ j with j an 
integer, and 
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where G\ — 2ir/a, with a the period of the UCDW. In- 
serting this into Eq. (0) we find E ucdw (Gi), the average 
energy per electron in the UCDW (see Fig. [j]). The op- 
timal UCDW corresponds to the minimum E ucdw , and 
is observed at a ~ 2.84Zfc\/2L + l (in general agreement 
with Ref. |7J], even though we are far from L— >oo), where 
each electron gains one to a few degrees, see Table Q. 
Since the anisotropic-isotropic transition is observed at 
temperatures much smaller than this, it is clear that the 
observed transition is not related to the formation of the 
stripes but, as we shall see, to the unbinding of topolog- 
ical defects in the stripes. 

(II) Low energy excitations of the UCDW. — Here we 
consider low energy states which correspond to long 
wavelength fluctuations of the UCDW. We take care to 
construct modulations of the stripes which do not accu- 
mulate charge over large distances since this would sig- 
nificantly increase the Coulomb energy of the system. 
These modulations add extra "Bragg peaks" to the den- 
sity function A(G) (see Fig. 0), and of the many modu- 
lations one can devise, very few avoid adding significant 
peaks far from where U(G) is near its minimum p2[ . 
These can be described by a distortion in the position of 
the UCDW stripe edges of the form 
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FIG. 2. Two characteristic examples of low energy per- 
turbations of the UCDW. Top: the longitudinal modulation. 
Bottom: the transverse modulation. On each panel, the 
right-hand side shows the Bragg peaks of A(G) in recipro- 
cal space. Gi is the wavevector of the UCDW, and q x , q y are 
the wavevectors of the modulation. See Eqs. (H)— (lid). 



u(x,y) = a cos(q x x) cos(q y y), 



(8) 



where a, q x , q y are the amplitude and wavevector com- 
ponents of the modulation respectively. Longitudinal 
(q y = 0) and transverse (q x = 0) modulations are illus- 
trated in Fig. 0. To determine the energy of this excited 
state to 0[a 2 ], we need to retain the following peaks: 
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where j is an integer. The energy per electron, relative 
to the optimal UCDW is then given by 
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Keeping terms up to 0[q x ,q y ,q x q 2 ], the energy per unit 
area is 

*£ = Y [Bql + Kq y + K'q 2 x q 2 y + K" q x ] , (12) 

with the elastic coefficients given by 
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TABLE I. UCDW: optimal wavevector Gi, period a, en- 
ergy gain per electron E ucd '" and elastic constants B and K. 
The calculations were performed for the realization of Ref. Ilf . 
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It is easy to see from energetics above [Eq. (O)] that 
the low-energy perturbations of a UCDW correspond 
one-to-one to those of a smectic liquid crystal [H] : 

E sm = lJd 2 r{[B {d x uf + K (d 2 y u) 2 ] 

+ [K'{d x d y uf + K"(dluf]}. (17) 

Results for the elastic moduli B and K are presented in 
Table Q for parameters relevant to the sample used in Ref 
IJJ. The terms between the second sets of brackets in 
Eq. (|1J]) are not expected to be relevant since they only 
become large for momenta near the edge of the Brillouin 
zone (where the validity of the elastic theory is doubtful) . 

We now use the energy functional E srn (without the 
terms involving K' and K") for all further analysis of 
the quantum Hall liquid crystal. 

(Ill) Effects of thermal fluctuations: from smectics to 
nematics. — The energy functional for a smectic [Eq. (17)] 
has been extensively studied. We follow closely the for- 
mulation of Toner and Nelson H|. Since the dimen- 
sionality of the system [d = 2) is one below the lower 
critical dimension for layered materials, phonon fluctu- 
ations readily destroy positional order for T > (the 
Landau-Pcicrls argument), while preserving order in the 
layer orientation. However, this argument omits dislo- 
cations, which have finite energy; their energy can be 
estimated as |22| 
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where A 2 = K/B and q c ~7r/a is a large momentum cut- 
off. Therefore, for T > we expect a density of disloca- 



tions given by nn 



than £d 
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, and as long as Ejj ^> fc^T, dislocations 



can be treated in a Debye-Huckel approximation. Then, 
to lowest order in q 2 and q 2 the correlation function for 
the layer normal angle 9 = —d y u can be written as |l4| 
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which is precisely the correlation function of a two- 
dimensional nematic, with a free energy 
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where n = (cos 9, sin 0) is the director field, and the two 
Frank constants are given by 



K x = K and K 3 = 1E L 



(21) 



Orientational correlations in the director n(r) should de- 
cay algebraically at distances much larger than £d ■ Table 
IJ summarizes the values of K\ and K%. The values of 
these elastic constants are determined at distances com- 
parable to £_ D (~ 10 a at T - 100 mK). 

(IV) The nematic to isotropic transition. — At suffi- 
ciently long wavelengths Nelson and Pelcovits pa], us- 
ing a momentum-shell renormalization approach, have 
shown that deviations from the one-Frank constant ap- 
proximations K\ — K3 are irrelevant, and the system is 
equivalent to a two-dimensional XY model: 



F xy = ±K(T)Jd 2 r(V6) 2 



(22) 



with K — > [i£i(£o) + K3(£d)]/2 at very large distances. 
For our values of K\ and K3, at the characteristic tem- 
peratures of the experiments, convergence is achieved at 
distances around 20-100 £d. We then expect unbinding 
of disclination pairs at the KT temperature |la] : 



IsbTkt = -3- K(Tkt) 



(23) 



where the tt/8 comes instead of the more common 7r/2 for 
vortices since each disclination winds up the angle by it 
rather than 27r. In general, K{Tkt) corresponds to the 
large-distance elastic constant (reduced by disclination 
pairs) to the bare elastic constant at small distances -ftT(O) 
by means of the KT RG formulas [15| : 
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where k = K/ksT and we have introduced the fugacity 
y ~ exp[— 7T 2 K(0)/kBT]. In practice, these RG equa- 



lly) tions can be approximated by ksTxT — (tt/8) K (0)/(l- 



27rexp[-7r 2 J fs:(0)/8fc B TKT]) ^ 0.86 (tt/8) if(0). This re- 
duction is in general agreement (although somewhat less 
important) to results for Monte-Carlo simulations |12] ]. 

Table O presents the resulting estimates for the discli- 
nation unbinding transition temperatures for half-filled 
LLs. Although these can only be considered estimates 
due to the approximations used, they are in quali- 
tative agreement with the temperatures at which the 



TABLE II. Frank elastic constants K\ and K3, renormal- 
ized elastic constant K and KT disclination unbinding tem- 
perature calculated for the experimental realization of Ref. 
[p|. Note the characteristic oscillations with the spin index. 
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anisotropies are seen to vanish. We also see the character- 
istic spin oscillation of the transition parameters |l|,p4|. 
The reason for this spin oscillation is simple: in the en- 
ergetics of Eqs. (|]-||), there is an energy scale e 2 /h that 
decreases with increasing filling factor v; simultaneously 
the matrix elements of the Coulomb interaction [Eq. (0)] 
increase with increasing LL index L, resulting in the ob- 
served spin dependence. 

There are a couple of caveats which apply to our re- 
sults. First, we've left out the native anisotropy of the 
sample which tends to align the smectic structure (similar 
effects arise from an in-plane component of the magnetic 
field). Uniaxial anisotropy will produce a term of the 
form B'(d y u) 2 in the smectic energy density; although 
the experiments indicate that B' <C B, at sufficiently 
long length scales (of order -J ' K/ 'B') the anisotropy will 
dominate over the bending energy. In this case the dislo- 
cation energy diverges as the logarithm of the system size, 
and the transition to the isotropic phase occurs through 
the unbinding of dislocations. Second, as is customary in 
studies of smectics, we have dropped terms in the smectic 
free energy, Eq. (|17|), of 0[q^q 2 ]. To check the validity of 
this truncation we have calculated the elastic coefficients 
K' and K" , and find that while K" > 0, it is possible for 
K' to be negative Eq] . This does not seem to cause any 
problems in the long wavelength limit, but it may change 
our estimates of the dislocation energy. This issue is cur- 
rently under study E3] . 

In conclusion, we have mapped a 2DES with half- 
filled LLs to a liquid crystal with smectic/nematic or- 
der at short/long distances and which undergoes a KT 
disclination unbinding transition, after which the system 
becomes isotropic, as seen by transport measurements. 
Without the use of any fitting parameters we have ob- 
tained transition temperatures in qualitative agreement 
with experimental evidence. A particularly robust fea- 
ture is the spin dependence of the nematic elastic moduli 
and transition temperature (Table O): they are larger for 
the lower spin sub-band (y = 9/2, 13/2, 17/2). While 
precise experimental values for the transition temper- 
atures have not been established and the transition is 
rounded by disorder, the same characteristic spin depen- 
dence is observed in the transport anisotropy p X x/pyy 
fiJig (see also Ref. @). 
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